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The objective of this paper is to present a simple pyrolysis model to capture the main characteristics 
of the decomposition of a thermally thin particle at high temperatures corresponding to those found in 
the furnace of coal/biomass power plants. To achieve this goal, it is assumed that pyrolysis begins soon 
after the surface of the particle has reached a certain pyrolysis temperature, and proceeds according 
to a shrinking (unreacted) core model with an infinitesimal reaction front. The formulation of various 
stages including initial heating, pre-pyrolysis heating, pyrolysis and post-pyrolysis heating is carried 
out based on a time and space integral method which allows one to describe the energy conservation 
equation in an algebraic form. Two different treatments are presented for the pyrolysis stage. The first 
formulation assumes separate temperature profiles for char and biomass regions (double-temperature 
profile), whereas in the second treatment only one profile is considered for the temperature throughout 
the particle (single-temperature profile). Of particular interest is the latter approach that leads to simple 
relationships for predicting the duration of various stages, enabling one to predict the mass loss history. 
The accuracy of both methods is examined by comparing their predictions with recent experimental 
data reported in the literature as well as the prediction of comprehensive pyrolysis models. Satisfactory 
agreement is achieved indicating that both pyrolysis models based on double- and single-temperature 
profiles can be used with sufficient accuracy for engineering purposes. 

© 2012 Elsevier B.V. All rights reserved. 


1. Introduction 

Thermochemical conversion of biomass is a topic of technolog¬ 
ical importance due to the increasing interest in the utilization 
of biomass as a renewable energy source. The main mechanism 
for chemically converting the biomass into other byproducts is 
heat from an external source. The energy from biomass can either 
be directly converted into heat, for instance, in the furnace of a 
power plant, or into the liquid fuel and biogas within a gasifier. 
The well-established methods of thermochemical conversion of 
biomass for engineering applications include combustion, gasifica¬ 
tion, torrefaction and pyrolysis. The key differences between these 
methods are related to the oxygen concentration in the surrounding 
fluid and the operating temperature. 

A substantial amount of research has been conducted on the 
subject of biomass pyrolysis by means of experimental and numer¬ 
ical investigations to increase the knowledge of various physical 
and chemical processes involved in the biomass pyrolysis. The 
reader may access the recent state-of-the art in Refs. [1-12]. The 
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common feature of the models reported in these references is that 
the transport equations are usually described one-dimensionally 
and combined with a kinetic scheme representing chemistry of the 
virgin biomass decomposition. Two and three-dimensional models 
have rarely been reported; see e.g. Ghabi et al. [13]. The advan¬ 
tage of these models is that they allow calculating a variety of 
parameters such as solid and gas phase densities, temperature, 
pressure and velocity fields within the porous matrix and forma¬ 
tion of main byproducts (light gases, tar and char) in space and time. 
Such advanced models enable one to get a deeper understanding 
of the physics and chemistry of the pyrolysis process on the scale 
of a single particle. 

Because of the transient nature of the pyrolysis process, the 
transport equations are necessarily described in partial differen¬ 
tial forms in order to capture the variation of key parameters along 
the particle and during the process. This would be at the expense of 
considerable amount of computational and programming efforts, in 
particular, when dealing with design of combustors and gasifiers, 
where a large number of particles undergo thermochemical conver¬ 
sion. For practical applications, designers are commonly interested 
in only a few parameters such as average temperature and/or 
temperature at the surface of the particle, rate and the amount 
of particle mass loss, ignition time at which a particle begins to 
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Nomenclature 

c p specific heat (J/kg K) 

h convective heat transfer coefficient (W/m 2 1<) 

k thermal conductivity (W/m I<) 

L thickness/Length of particle (m) 

L cr minimum particle length for transition from ther¬ 
mally thin to thermally thick regime (m) 
m" decomposition rate per unit surface area (kg/m 2 s) 
q" heat flux (W/m 2 ) 

q" et net heat flux at particle surface (W/m 2 ) 

R particle radius (m) 

T temperature (I<) 

T p pyrolysis temperature (K) 

T s surface temperature (K) 

T sL surface temperature at the time t L (K) 

t time (s) 

ti duration of thermal penetration movement from 

the front surface to the back face of the particle (s) 
tp i n i time of commencement of pyrolysis at the surface 
of particle (s) 

t p conversion/pyrolysis time (s) 

x spatial coordinate 

x c char depth (m) 

x t thermal penetration depth (m) 

Greek letters 

a thermal diffusivity (m 2 /s) 

Ah p specific enthalpy of pyrolysis (J/kg) 

s emissivity 

p density (kg/m 3 ) 

a Stephen-Boltzmann coefficient (W/m 2 I< 4 ) 

tp coefficient of Eqs. (3) and (17) 

coefficients of Eq. (16) 

Subscripts 

0 initial condition 

oo surrounding condition 

B biomass 

BC boundary condition 

C char 

ext external 

f final 

ph post-pyrolysis heating up 

r reactor 


decompose, and total duration of the conversion process. To cap¬ 
ture the main characteristics of the pyrolysis of solid fuels with 
optimized computational efforts requires one to employ less com¬ 
plex models. 

From a mathematical point of view, a key step towards simpli¬ 
fied pyrolysis models is to describe the particle conversion process 
in the form of a set of ordinary differential equations (ODEs) rather 
than partial differential equations (PDEs). In principle, a chemically 
converting solid particle may ultimately decompose in the follow¬ 
ing two extreme limits: the regime of shrinking density, and the 
shrinking core regime. In the former limit, a particle undergoes an 
almost homogenous conversion in the absence of intra-particle gra¬ 
dients, whereas in the latter extreme limit the reaction takes place 
at a very thin layer, resulting in a reaction front that is created at 
the surface of the particle and moves towards its center. 

For the case of pyrolysis of charring materials, e.g. wood, a 
particle may decompose homogeneously in the regime of kinet- 
ically controlled conversion. Our previous study [2] indicated 


Temperature [K] 


1200 



0.015 


Char density [g/cm3] 



0.005 


Fig. 1 . Illustrative time and space evolution of intra-particle gradients of a biomass 
particle undergoing a homogenous pyrolysis (250 |jim; 1050 K; 300 kg/m 3 ). Different 
lines correspond to different times increasing in the direction of the arrows. 


that homogenous pyrolysis of a wood particle may occur at low 
densities, small particles and low operational temperatures corre¬ 
sponding to (typically for a cylindrical beech wood) Bi<0.1 and 
Py’>1.5 x 10 -3 , where Bi and Py’ denote Biot and the external 
pyrolysis numbers, respectively. A typical time and space evolution 
of biomass and char densities as well as temperature of a particle 
undergoing a homogenous pyrolysis predicted by a comprehensive 
pyrolysis model [1,2] is depicted in Fig. 1. 

It is however unlikely that pyrolysis of biomass in combustors 
and gasifiers occurs in the regime of kinetically controlled in the 
absence of intra-particle gradients, due to the high operational tem¬ 
peratures and relatively large particles. It is difficult to mill biomass 
due to its fibrous structure which results in particles of relatively 
large size (around 1 mm). Moreover, Lu et al. [14] and Bharadwaj 
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et al. [15] examined the accuracy of a lumped model (treating the 
particle as a whole) and found significant errors in the predictions 
indicating that the intra-particle effects need to be accounted for 
in a pyrolysis model. 

Based on the arguments provided in Refs. [2,14,15], one needs to 
account for the temperature gradients inside a pyrolyzing particle. 
A recent study by Saastamoinen et al. [16] has also revealed that 
the intra-particle gradients of particles of less than a millimeter 
in real furnaces can be significant. The pyrolysis model employed 
by Saastamoinen et al. [16] assumes that the density of particle 
decreases homogeneously throughout the particle. To account for 
internal temperature gradients their model incorporates an inter¬ 
nal thermal resistance in terms of particle thermal conductivity and 
size. 

Another case is to employ the concept of shrinking core model 
assuming that the pyrolysis begins as soon as the surface tempera¬ 
ture reaches a characteristics temperature T p (to be called pyrolysis 
temperature throughout this article), thereby yielding a reaction 
front. A further simplifying assumption is that the virgin biomass 
decomposes to char and volatiles at an infinite rate in a very thin 
layer dividing the particle into char layer and biomass region. To 
account for the temperature gradient inside the particle, one may 
employ the concept of integral method that was originally pro¬ 
posed by Goodman [17]. This method allows one to describe a 
one-dimensional transient heat transfer problem represented in 
terms of a partial differential equation (PDE), in the form of an 
ordinary differential equation (ODE). 

The integral method is based on an assumption for the spatial 
behavior of the temperature along the particle, so that upon inte¬ 
grating the energy equation over the space coordinate the particle 
pyrolysis model is converted into an ODE. This method has been 
used for predicting the pyrolysis of solid particles in some past 
studies [18-26]. Although the concept of the shrinking core model 
has been commonly employed in these studies, the essential dif¬ 
ferences between them are summarized in Table 1. Furthermore, 
the formulation and solution methodology differ from one study to 
another. For instance, an infinite decomposition rate at a constant 
pyrolysis temperature assumption is adopted in Refs. [18-24]. On 
the other hand, Galgano and Di Blasi [25,26] used a finite rate model 
of conversion for cylindrical wood particles. 

Given that the concept of a fixed pyrolysis temperature T p may 
seem to be controversial, it is a useful simplifying approach, in par¬ 
ticular, when dealing with flaming combustion [27]. Furthermore, 
Galgano and Di Blasi [28] assessed the predictions of a shrinking 
core model with both finite rate and infinite rate pyrolysis models. 
They found that the unknown temperature T p is not necessarily 
the same as that predicted by the finite rate model. However, the 
main conclusion was that the most important process characteris¬ 
tics could be accurately predicted by assigning a proper value for 

V 

A careful review of the above cited sources indicates that the 
concept of shrinking core model has been only used for predict¬ 
ing the pyrolysis of thermally thick particles in which ignition takes 
place at the particle surface before thermal penetration reaches 
the back face. A further possible situation which has not been fully 
realized in past studies is the case in which the ignition occurs after 
thermal penetration has reached the back face. In this study, this 
will be referred to as thermally thin particle. More importantly, the 
various stages of conversion of a finite size particle have not been 
adequately described. These observations form the primary moti¬ 
vation of the present study to develop simple models for a finite size 
thermally thin wood particle undergoing a pyrolysis process. The 
main objective is to introduce a novel modeling approach based on 
a time and space integral method , which can be easily linked with 
CFD models for designing combustors and gasifiers of solid par¬ 
ticles. It will be shown that the model presented in this paper is 






Fig. 2. Schematic representation of a pyrolyzing thermally thin particle. 


much simpler and easier for coding compared to those presented 
by others. 

In the upcoming section, we first explain the sub-processes 
involved in the pyrolysis of a thermally thin particle. The formu¬ 
lation of various stages will be subsequently presented in Section 
3. Typical numerical results and the validation of the new model 
against different experiments and the prediction of a comprehen¬ 
sive model will be given in Section 4. Finally, Section 5 deals with 
the evaluation of the model results and its accuracy. 

2. Description of the process 

Consider a dry slab biomass particle with a thickness of L that 
is exposed to a hot environment. The particle is initially at tem¬ 
perature T 0 (<T P ), and it is assumed that the pyrolysis begins as 
soon as the surface temperature reaches the pyrolysis tempera¬ 
ture T p . Further assumptions taken into account in the development 
of the present model (Section 3) are: (1) particle volume remains 
unchanged during the entire process; (2) thermo-physical prop¬ 
erties are constant; (3) pyrolysis takes place at a very thin layer 
where the temperature is T p ; (4) the convective flow of volatiles is 
neglected. 

Depending on particle size, thermo-physical properties and 
external heating rate, the particle may experience two different 
conversion regimes as depicted in Figs. 2 and 3. Initially, the par¬ 
ticle begins to heat up so that a thermal wave with a thickness x t 
measured from the surface is formed. One possible situation is that 
the thermal wave reaches the back face (assumed to be insulated in 
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Table 1 

The specification of the integral method employed in past studies. 
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Researched s) 

Temperature profile 

Particle 

Validation against 

Kanury [18] 

Linear 

Charring slab and cylinder 

Experiments 

Chen et al. [19] 

Exponential 

Charring and non-charring finite slab 

Analytical solution 

Quintiere and Iqbal [20] 

Quadratic 

Non-charring semi-infinite slab 

Experiments and PDE model 

Weng et al. [21] 

Linear (char layer) and quadratic (virgin layer) 

Charring semi-infinite slab 

Experiments of wood pyrolysis 
with char oxidation 

Spearpoint and Quintiere [22] 

Linear (char layer) and quadratic (virgin layer) 

Charring semi-infinite slab 

Experiments of different 
pyrolyzing woods 

Galgano and Di Blasi [24] 

Quadratic 

Charring cylinder 

Experiments of wood pyrolysis 


this study); i.e. x t equals L, while the surface temperature T s is still 
less than T p (see Fig. 2b). The heating up process continues until T s 
becomes T p ; at this moment the surface of the particle decomposes 
to volatiles and char residue. On the other hand, a second possibil¬ 
ity is that the thermal wave is still moving towards the back face of 
the particle while T s has already reached T p (see Fig. 3b). In order 
to avoid any confusion, we intentionally distinguish between these 
two situations; the first case will be referred to as thermally thin 
particle , and the second one as thermally thick particle. 

A thermally thin particle undergoes the following four stages: 


1. Initial heating up (Fig. 2a), 

2. Pre-pyrolysis heating up (Fig. 2b), 

3. Pyrolysis (Fig. 2c), 

4. Post-pyrolysis heating up (Fig. 2d). 






Fig. 3. Schematic representation of a pyrolyzing thermally thick particle. 


In the first stage of the process, a thermal wave (Fig. 2a) is cre¬ 
ated at the surface and it moves inside of the particle as the time 
passes until it reaches the back face; i.e. x t = L, at the time t L while 
the surface temperature is still less than T p . As shown in Fig. 2b, 
the heating up process continues until T s reaches T p ; whereby ini¬ 
tiating the pyrolysis process at the surface of the particle. The time 
required for the surface of the particle to reach the pyrolysis tem¬ 
perature is denoted by t p ini . From this moment on, the particle is 
divided into two regions: the virgin part, which still undergoes the 
heating up process, and the char layer with thickness x c (Fig. 2c). The 
pyrolysis front propagates inside the particle by continuation of the 
process so that the thickness of the char region increases until time 
t p at which x c equals L. At this moment the virgin biomass has com¬ 
pletely converted to char and volatiles assumed to leave the particle 
immediately upon formation. The final stage is the char heating up 
(post-pyrolysis stage) as depicted in Fig. 2d. So, the problem is a 
simple conduction heat transfer. This stage will cease when thermal 
equilibrium between the particle and surrounding is established. 

As pointed out earlier, for a thermally thick particle the pro¬ 
cess is a little different (see Fig. 3) in a way that after initiation of 
the heating up stage (Fig. 3a), the surface temperature increases 
up to T p before the thermal wave reaches the back face. In other 
words, for a thermally thick particle ti > t pini . Thus, upon initiation 
of the pyrolysis process, the particle is divided into three regions: 
char layer (0 < x < x c ), virgin biomass undergoing heating up process 
(x c <x<xt), and virgin biomass maintained at the initial tempera¬ 
ture (x t < x < L). By continuation of the pyrolysis process, the thermal 
wave and the char front penetrate into the particle until x t reaches 
the back face. Beyond this moment, the rest of the process is the 
same as the thermally thin particle as depicted in Fig. 3c and d. In 
summary, the following processes take place in a thermally thick 
particle: 

1. Initial heating up (Fig. 3a), 

2. Fleating-pyrolysis (Fig. 3b), 

3. Pyrolysis (Fig. 3c), 

4. Post-pyrolysis heating up (Fig. 3d). 

The focus of the present study is on thermally thin particles. 
In the upcoming section, the formulation of the four stages of the 
pyrolysis process will be presented using a time and space integral 
method. 

3. Formulation 

The main basis for model development is that the temperature 
within the particle will be assumed to obey a quadratic function sat¬ 
isfying the boundary conditions. The assumption of decomposition 
of virgin biomass in an infinitesimal thin layer whose tempera¬ 
ture is represented by T p leads to simply describe the heat transfer 
equation in terms of transient conduction problem in both char 
and biomass layers. Upon integrating the heat transfer equation 
with respect to the space coordinate, the PDE form of energy equa¬ 
tion reduces to an ODE with respect to time. In the next step, time 
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integration is applied to the resulted ODE which leads to an alge¬ 
braic equation. This approach will be used in all stages of the 
conversion process depicted in Fig. 2. 

3.1. Initial heating up 


As the particle undergoes a transient conduction process 
(Fig. 2a), the conservation of energy can be represented as follows: 


dT 


d 2 T 


dt ttB dx 2 ' 


(1) 


where a B = /<b/(PbC p b). 

The following initial and boundary conditions are prescribed for 


this problem: 



Initial condition : T 

= T 0 at 

[ = 0 . 

Boundary condition : 

EC, : 

k 9T i 

^B ^ x=0 — <7 net 


BC 2 : 

k B §U=* t =°, 


BC 3 : 

T\ x=xt = T 0 . 


For a known external radiant heat flux q" xt , the net heat flux at the 
surface of the particle is: 


<?net = Qext ~ h ( T s ~ Too) - cre(T 4 - T^). 


( 2 a) 


On the other hand, for a known reactor temperature T r , the net heat 
flux is described as: 


q'net = 'l(Tr-T s ) + ( Te(T r 4 -T s 4 ). 


( 2 b) 


As mentioned earlier, the temperature profile is assumed to obey a 
quadratic function at each instant. Hence, 


T — (p 2 X 2 + <p\X + (po, 


(3) 


where the coefficients of Eq. (3) are functions of time. With the aid 
of the boundary conditions mentioned above, these coefficients can 
be obtained so that the resulting temperature profile obeys: 




(4) 


Applying a space integration to Eq. (1) from x = 0 to x = x t yields: 


d_ 

dt 


Tdx-T 0 -^ ) -a B [ — |x=x t - ^ 7 lx=o 


dT 

a^ 1 


dT 

a^ 1 


(5) 


( 6 ) 


Using Eq. (4) and the boundary conditions, Eq. (5) reduces to: 

^(Qnet^t ) = ^Ot B q net 

Applying approximate time integration to Eq. ( 6 ) using a mean 
value for q" et (= [q" et + q"] / 2 ) leads to: 


x t (t)= 4 / 3Q!b 


1 + 


Ct 


(7) 


where q" denotes the net heat flux at t = 0. Note that the net heat flux 
is a function of T s which is obtained from the temperature profile 
relationship; i.e. Eq. (4). Hence, 


r, (t ).r 0 + 2 g„. 


( 8 ) 


The time required for the thermal wave to reach the back face 
is obtained by inserting x t = L into Eq. (7). Hence, 


ti = 


L 2 


3«b[ 1 +(<7o/<et)l 


0) 


For a thermally thin particle, t L is less than t p ini . At a given pro¬ 
cess condition, it is possible to determine the particle thickness at 
which t L equals t pini . To determine the time at which T s reaches T p ; 
i.e. initiation of pyrolysis, one should first eliminate x t between Eqs. 
(7) and ( 8 ), rearrange the resulting expression for t and replace T s 
by T p . Hence, 


(Tp ~ 7p) 2 _ 

3a B (q" et /2k B f^+(q''/q" et ))' 


( 10 ) 


Equating the right hand sides of Eqs. (9) and (10) and rearranging 
for L gives: 


Ter — 2 /<b 


( Tp-T 0 \ 

\ ^net / 


( 11 ) 


L cr denotes the minimum size of the particle for transition from 
thermally thin to thermally thick regime. Notice that for a particle 
heated uniformly in a reactor from both sides, L (and L cr ) would 
exhibit half of the particle thickness. 

In fact, Eq. (11) allows one to identify whether a particle is ther¬ 
mally thin or thermally thick. If the particle size is less than L cr , the 
problem under study is characterized with the regime of thermally 
thin otherwise the particle is thermally thick. It should be noted that 
for a thermally thin particle (L<L Cr ), t Piini has to be obtained from 
the next stage of the formulation to be followed in the subsequent 
section. 

Typically, thermal conductivity of various biomass types is 
between 0.2 and 0.4W/mK, the pyrolysis temperature varies in 
the range of 523-773 K, and the reactor temperature is above 
1200 K. For instance, assume T p = 673 K, T 0 = 300 K, k B = 0.25 W/m K, 
T r = 1273 K, h = 15 W/m 2 K. Using Eq. (2b), the net heat flux at the 
surface of the particle at the time of initiation of pyrolysis is 146 kW. 
Inserting the above values in Eq. (11) gives L cr = 1.3 mm, which is 
equivalent to a particle thickness of 2.6 mm if the particle is heated 
from both sides. In practice, such as in power plant furnaces, the 
particle thickness is usually less than 1 mm for raw biomass parti¬ 
cles. In the case of torrefied biomass, the particle size is expected 
to be even smaller; e.g. less than 0.5 mm. This sample calculation 
indicates that the regime of particle pyrolysis in real combustors 
is thermally thin (Fig. 2) since the particle thickness is less than its 
critical value. It is obvious that depending on ignition temperature, 
thermal conductivity and heating condition, the value of L cr would 
be different. Typical values of the critical half thickness of a particle 
are depicted in Fig. 4, which supports the idea that the pyrolysis of 
(raw and torrefied) biomass particles in most real furnaces occurs 
in a thermally thin regime. 


3.2. Pre-pyrolysis heating up 

In the second phase of heating up process, the thermal wave 
has already reached the back face while T s is still less than T p (see 
Fig. 2b). A similar procedure as explained above is undertaken. That 
is, the problem is still described by Eq. (1) assuming that the tem¬ 
perature profile can be approximated with Eq. (3). The following 
initial and boundary conditions are applied: 

Initial condition : T s = T s L at t = t L . 


Boundary condition : BCi : - /c B |*=o = q'nev 

BC 2 • Tb-^-Ix=l = 0, 

BC 3 : T | x=0 = T s . 
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T P [K] 



Tp(K] 


Fig. 4. Critical half thickness of particle versus ignition temperature at varying ther¬ 
mal conductivity at a reactor temperature of (a) 1273 K, and (b) 1473 K. 


Upon determining the coefficients of Eq. (3) with the aid of the 
above boundary conditions, the temperature takes the following 
form as a function of x: 


T = T S — 


l< B 


^net y 2 

2k B L 


( 12 ) 


Integrating Eq. (1) with respect to x along the particle thickness, 
and inserting Eq. (12) into the resulting expression yields: 


To determine the time at which the pyrolysis initiates at the 
surface, one needs to solve Eq. (14) for t with T s = T p . This leads to: 


fp.ini — f L + 


T P - T sL - (q" ini - q") — 


2 k B L 


3h\ a B (q;. ni + q'()’ 


(15) 


where q" ini is the net surface heat flux at the time of initiation of 
pyrolysis to be calculated using Eq (2) with T s = T p , T sL denotes the 
surface temperature at ti which can be determined from Eqs. (2) 
and (8) with x t = L, and qis the net surface heat flux at the time t L 
to be obtained from Eq (2) with T s = T sL . 


3.3. Pyrolysis 

The formulation presented in the previous section is valid until 
initiation of pyrolysis at the surface of the particle. According to 
Fig. 2c, the particle consists of a char region with thickness x c and a 
virgin biomass region with length L-x c . The temperature at the vir¬ 
gin/char interface, where the pyrolysis takes place at infinitesimal 
rate, is T p . Consistent with the analysis presented for the heating 
up stage, the temperature inside the particle is assumed to obey 
a quadratic function. Two treatments will be presented: Double¬ 
temperature profile, and single-temperature profile. In the former 
case, separate temperature profiles are considered for biomass and 
char regions, whereas in the case of single temperature profile, 
only one spatial profile is assumed for temperature throughout the 
particle. 


3.3.1. Double-temperature profile 

Let us first consider the case of double-temperature profile such 
that the temperature profiles in the char and biomass layers are 
represented as follows: 

Char region(0 < x < x c ): T c = ir 2 {x c - x) 2 + 0i(x c - x) + 


(16) 


Biomass region(x c < x < L): T B = 02(* - *c) + 0i(x - x c ) + 0o- 

(17) 

The heat transfer equation in the biomass region is still represented 
by Eq. (1). A similar equation can be prescribed for the char layer. 


8T 8 2 T 

9 F - “ c 9 x 2 ’ 
where a c = fec/(/OcC p c). 

By applying the following boundary conditions, 
BC. : -kc|^l*=o = <?net> 


(18) 


Rr • 9Tb i n 

bc 2 .^i x=l = 0, 


BC 3 : TbIx=x c — Tclx=x c — Tp, 


d_ 

di 



^net j 2 

3/<b 



(13) 


i 9Tg 8T(2 ■ rr a 1 

BC 4 : -/<B-^lx=x c = -kc-^-\x=xc +m"Ahp, 


Integrating Eq. (13) with respect to the time from t L to t, and 
assuming a mean value for q" et equaling (q" et + q")/2 ) leads to the 
following relationship: 

Ts = T sX + (q'net - + 2 _ tl ^ 


and determining the coefficients i/q and 0/, Eqs. (16) and (17) are 
rewritten as: 


Tc = T p - 


k B (p 1 + rh"Ah p q'^ et + k B </> 1 + m" Ah p ^ 


k c 


- (x c — x) + 


2x c k c 


(x c — X), 
(19) 


where the surface temperature and the net heat flux at the time t B , T B = T + 0 t (x - x c ) - ^ ( x _ Xc ) 2 

i.e. T s L and q'^ are obtained from the previous heating up phase. 2 (L - x c ) 


( 20 ) 
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where Ah p denotes the specific heat of pyrolysis (negative when 
endothermic and positive when exothermic), and m" represents the 
decomposition rate per unit particle surface area perpendicular to 
x-coordinate. 


Rearranging the above equation and performing approximate time 
integration between t pini and t, we get: 

PcCpCQ net(^^c — 2x c ) — 2p^Cp^q ne ^{L — ^PB^pBQp jni^ 


m" = p B 


Ax c 
AF’ 


( 21 ) 


= 3fccWnet + + Qp ini )(t - t p ,ini). 


(30) 


where At is a very small increment of time. 

Initially; i.e. t = t p ini ,x c = Oand0i = -q"//^. Integrating Eq. (18) 
with respect to x between x = 0 and x = x c , and using boundary con¬ 
ditions and temperature profiles given in Eq. (19) yields: 

^[(<Jnet - 2fc B0i - 2m"Ahp)x 2 C ] = 6ac(q"et + Mi + m"Ah p ). (22) 
Approximate time integration of Eq. (22) from t p ini to t gives: 


x c 



( Ct + Mi +m"Afi p V 

— 2/c b </>i - 2th" Ahp J p ' ini 


(23) 


Note that a mean value for the terms inside the bracket on the 
right hand side of Eq. (22) is considered when performing the time 
integration. 

Likewise, integrating Eq. (1) between x = x c and x = L using Eq. 
(20) results in: 

^ [01 (L - Xc) 2 ] = -3 « b 0i - (24) 

After applying approximate time integration from t- At to t and 
some rearrangements, a solution of Eq. (24) over a small time incre¬ 
ment At is found: 


0i(O = 0i(t- At) 


L - x c (t - At)' 

2 

3q:bA t 

. L-Xc(t) . 

exp 

. ( l-Xcf. 


(25) 


where x c = [x c (t) + x c (t - At)]/2. 

At each instant, the surface temperature is determined from Eq. 
(19) withx = 0. Hence, 


Ts=Tn + 


CfMi - m"Ah p 
2 k c 


x c . 


(26) 


3.3.2. Single-temperature profile 

Similar to the analysis presented in Sections 3.1 and 3.2, we 
assume that the temperature profile throughout the particle is rep¬ 
resented by Eq. (3) whose coefficients can be obtained with the 
appropriate boundary conditions at the particle surface and back 
face, and at the location of x c where the temperature is T p , thereby 
resulting in: 

r - r --% £t T £(x - ,i ‘> + < 27 > 

From conservation of energy for the particle, the net change 
in the energy of the particle is equal to the net amount of energy 
transferred to the particle. This would lead to: 

d [ l t , _ dx c d /%. dx c 

PBCpB-^j; J Tax + p B c P B7p-j^r + PcCpc^j. J Tax - p c c p cT p — 

= Qnet + th /y 'Ahp. (28) 


Substituting Eq. (27) into Eq. (28) and integrating yields: 


3 kcL dt [l? " et(i Xc)3] + 6 kcL dt 2x < 51 

= q" et + m"Ah p . 


( 29 ) 


Eq. (30) can further be reshaped to read: 


t = 

+ 


fp.ini 

PcCpcQnet( 3Lx c ~ 2 *D ~ 2 PBCpBq'nec( L ~ x cf + 2pBC p Bq" M L 3 


3k c L(q' I [ et + m"Ah p + q" i J 


(31) 


It is possible to estimate the pyrolysis time using Eq. (31). The 
pyrolysis process ceases when the char penetration depth x c equals 
I. At this moment, the mass loss rate is zero. Thus, inserting x c =L 
and m" = 0 into Eq. (31) leads to a simple relationship for estimating 
the pyrolysis time. 


£p — fp.ini + 


L 2 (.PcCpcQp + 2£> B c pB qp inj ) 
3k c(<7p + Qpjnj) 


(32) 


where q" denotes the net surface heat flux at the time of complete 
conversion of the particle. To determine its value, one needs to solve 
Eq. (2) together with an equation for the surface temperature at 
time t p to be obtained from Eq. (27) with x = 0 and x c = L, 


Q'v l 


(33) 


In Eq. (32), t pini can be well-estimated using Eq. (15), and q" 
is calculated from Eqs. (2) and (33). The rest of the parameters 
appearing in Eq. (32) are in fact treated as the problem input. In 
the forthcoming section, we will discuss that the usefulness of Eq. 
(32) is not limited to only approximating the pyrolysis time. In fact, 
it may further allow one to predict the mass loss history of a particle. 

Under special circumstances and depending on the operational 
conditions, the net surface heat flux may approach zero. In this case, 
Eq. (32) simply reduces to: 


, . 2l2 PbC p b 

L P ~ L P,ini + 3 / Cc 


(34) 


Eqs. (32) and (34) reveal that the pyrolysis time is proportional 
to the biomass density and to the square of the particle size. We 
arrived at a similar conclusion in a previous study [2] in which it 
was shown by means of a comprehensive pyrolysis model [1 ] that 
the total conversion time of a pyrolyzing biomass particle varies 
linearly with the initial particle density and it changes in a quadratic 
functional form with the particle size. From a physical point of view, 
the higher the heat capacitance of the biomass, i.e. PbC p b, the longer 
the duration of the char layer movement in order to completely 
penetrate into the particle. Indeed, a biomass particle with a higher 
heat capacitance requires longer time in order to absorb a certain 
amount of heat. On the other hand, it can be seen in Eqs. (32) and 
(34) that the pyrolysis time is proportional to the inverse of char 
thermal conductivity. This is because at higher values of kc, the heat 
received at the surface of the particle rapidly transfers through the 
char layer so the heat transfer rate into the biomass layer becomes 
faster. 


3.4. Post-pyrolysis heating up 

Upon completion of the pyrolysis, the particle becomes com¬ 
pletely char so that the problem reduces to a simple conduction 
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heat transfer problem (Fig. 2d). Therefore, the solution of the post¬ 
pyrolysis process would be identical to that of the second phase of 
biomass heating up as described in Section 3.2 with the difference 
that all thermo-physical properties should be replaced with those 
of char. Thus, the temperature profile and the surface temperature 
are obtained as follows: 


1351 

T s = Tsp + (q'net _ + “ f p)- (36) 

Eq. (36) can be rearranged as follows: 

r s -r sp -(q" et -q")(L/3kc) 

p (ac/2kcLXq''et + q£) ' 

At the final moment of the post-pyrolysis heating stage, the sur¬ 
face temperature is equal to reactor temperature T r , so there would 
exist thermal equilibrium between the particle and the surround¬ 
ing, implying a zero net heat flux at the surface of the particle. Thus, 
Eq. (37) would reduce to: 


fph — ff fp 


T r -T sp + (q''L/3k c L) 
(a c qp/2k c L) 


(38) 


where t ph denotes the duration of the post-pyrolysis heating up 
stage. 

Eliminating q " between Eqs. (33) and (38), we find: 




(39) 


It can be inferred from Eq. (39) that the minimum duration of 
the post-pyrolysis stage would occur if T sp = T r . Elence, 


fph — 


2 L 2 
3ac 


(40) 


It has been observed from detailed numerical models [2,4,15,16] 
that the trend of mass loss history of small biomass particles at high 
temperatures is almost linear. We will later show that the mass 
loss history predicted by the double-temperature profile model 
also approximately obeys a linear function. Since the final parti¬ 
cle density is treated as a known parameter in the present study 
(as well as in most simple pyrolysis models; e.g. a model based 
on the global one-step kinetics), one may locate the starting and 
ending points of mass loss curve in a mass loss-time plot, as illus¬ 
trated in Fig. 5. The initiation of the pyrolysis would correspond to 
a mass loss of zero and the time tp iini to be estimated using Eq. (15). 
Furthermore, the end point of the mass loss history graph would 
correspond to 1 - PcIpb, and the time t p to be obtained from Eq. 
(32). The mass loss versus time would simply be a line connecting 
the above two points on a mass loss-time plot. The last stage of 
the process; i.e. post-pyrolysis heating, is the location (1 -Pc/Pb» 
f P + tph) as depicted in Fig. 5. It will be shown in the next section 
that mass loss graphs obtained from this method are well compa¬ 
rable with those resulted from a comprehensive pyrolysis model. 
Given the simplicity of this method for predicting the mass loss his¬ 
tory, it does not provide a prediction of the mass loss rate history. 
Rather, it gives an average rate of mass loss flux throughout the pro¬ 
cess defined as (m" - m")/(t p - tpjni) where m" = pL represents the 
particle mass per surface area. Hence, 


f 3k c (p B -A:)(q;) + q',; iini ) 

ave UPcCpcQp + 2pBC P B<Jp ini ) ' 


(41) 


It is also possible to provide an approximate time evolution 
of the surface temperature with a similar procedure. As outlined 
above, the duration of various stages depicted in Fig. 2; i.e. t L , t p ini , 



Time 


Fig. 5. Schematic diagram for predicting the mass loss history based on the single¬ 
temperature profile approach. 


t p , tf can be determined using the appropriate relationships pre¬ 
sented in previous sections. On the other hand, the corresponding 
surface temperatures; i.e. T L , T p , T sp , T f , are either already known 
(T p and T f ) or can be determined. Thus, a line connecting (T 0 , 0), 
(T s l, t L ), (T p , tpi ni ), (Tsp, tp), (Tf, t f ) would provide an approximation 
of the surface temperature history (Fig. 6). 

3.5. Numerical solution 

To provide a convenient means for designers, the calculation 
procedure of a thermally thin particle using the proposed model 
with pyrolysis stage formulated based on the double-temperature 
profile approach is illustrated in Fig. 7. The advantage of this model 
compared to the one based on the single-temperature profile treat¬ 
ment is that it allows calculating the history of mass loss rate and 



Time 

Fig. 6. Schematic diagram for predicting the surface temperature history based on 
the single-temperature profile approach. 
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Fig. 7. Calculation flowchart for predicting the pyrolysis process based on the 
double-temperature profile approach. 


other key parameters by solving the appropriate equations of each 
stage at incremental time steps. In the initial heating up process, 
one needs to simultaneously solve Eqs. (2), (7) and (8). The calcu¬ 
lations continue until x t equals L at the characteristic time t L . Upon 
satisfaction of x t = L, Eqs. (2) and (14) are solved to compute the net 
heat flux and the temperature at the surface of the particle during 


the pre-pyrolysis heating up stage. When the surface temperature 
attains an assigned pyrolysis temperature at the time t p ini , the solu¬ 
tion algorithm switches so that Eqs. (2), (21), (23), (25) and (26) are 
solved numerically through a trial-and-error method until x c equals 
I. Finally, Eqs. (2) and (36) are solved until establishment of thermal 
equilibrium between particle and surrounding. 


4. Results and discussion 

The accuracy of the presented models based on an approximate 
time and space integral method has been examined by comparing 
the predictions of the model in terms of ignition/pyrolysis time; i.e. 
tp ini, mass loss history and the pyrolysis time against the experi¬ 
ments of different wood particles reported in the literature as well 
as the computations of a comprehensive model of pyrolyzing single 
wood particle published by the authors earlier [1,2], which was fur¬ 
ther used in a particle combustion model [29]. The thermo-physical 
properties employed in the simulations are those reported in the 
relevant sources taken for validation of the model, or given in the 
literature. Unlike in many past studies in which an endothermic 
heat of pyrolysis has been employed, in this study the correla¬ 
tion proposed by Milosavljevic et al. [30] is used for calculating 
the enthalpy of pyrolysis. It is proven to be a sufficient relation¬ 
ship, which allows producing accurate predictions when used in a 
pyrolysis model [2]. 

The first set of validation is related to prediction-experiment 
comparison of ignition time as given in Table 2. The data are 
obtained from Lu et al. [4] and Brescianini et al. [31 ]. In the experi¬ 
ments of Lu et al. [4], sawdust particles with equivalent thickness of 
320 |jim were heated up in an entrained-flow reactor maintained 
at 1625 K. Brescianini et al. [31] studied the ignition characteris¬ 
tics of plywood exposed to an incident heat flux in the range of 
25-50 kW/m 2 . It can be seen that the predicted values of ignition 
time compare well with the measured values. 

The measured mass loss history of the sawdust particles was also 
reported by Lu et al. [4], so these data are further used for valida¬ 
tion of the simplified pyrolysis models presented in this study. The 
comparison of the predicted mass loss histories using double- and 
single-temperature profile pyrolysis models with the measured 
values and the prediction of the comprehensive pyrolysis model 
of Lu et al. [4] is shown in Fig. 8. The results of both treatments 
compare very well with the model prediction of Lu et al. The trend 
of mass loss history and the conversion time reasonably agree with 
the experiments. All models exhibit a rather steeper mass loss his¬ 
tory compared to the experimental data, which is likely due to the 
imperfect size and irregularity in particle shape, as pointed out by 
Lu et al. [4]. 

In a previous study by Haseli et al. [2], the effects of particle size 
and external radiant heat flux on conversion time and final char 
yield of beech and spruce particles were studied by means of a com¬ 
prehensive pyrolysis model. The detailed description of this model 
can be found in Refs. [1,2]. The comprehensive model, validated 
against a wide range of experiments, allows one to observe time 
and space evolution of many parameters such as temperature, solid 
density, gas phase density, porosity, velocity of gaseous flow within 
the pores, etc. The conversion times of beech and spruce wood slab 
particles predicted by the simplified models of this study and the 
comprehensive model are illustrated in Figs. 9 and 10, respectively. 

Figs. 9a and 10a illustrate comparisons between these mod¬ 
els for an incident heat flux of 100kW/m 2 and varying particle 
size. Simulations were carried out for the situation that both front 
surface and back face of particle are exposed to the above heat 
flux, so that the calculations were performed for half of the par¬ 
ticle assuming a symmetry condition at the particle center. The 
comparisons shown in Figs. 9a and 10a reveal the capability of 
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Table 2 

Comparison of the predicted and measured ignition time. 
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Wood type 

Plywood 

Plywood 

Plywood 

Sawdust 


Size [mm] 

11.5 

11.5 

11.5 

0.16 


Density [kg/m 3 ] 

550 

550 

550 

650 


T P [K] 

673 

673 

673 

573 

External heating 

Q'ext [l<W/m 2 ] 

25 

35 

50 

- 


Tr [K] 

- 

- 

- 

1625 

tig [s] 

Measured 

113 ±20 

45 ±4 

22 ±6 

0.075 

Calculated 

133.8 

48.2 

20.1 

0.089 


Source of data 

Brescianini et al. (2003) 



Lu et al. (2010) 


the simplified models for predicting the conversion time of small 
particles. In order to indicate the effect of the external heating con¬ 
dition, another comparison is made in terms of the conversion time 
for a 1 mm particle at varying external heat flux, as depicted in 
Figs. 9b and 10b. As can be seen, the simplified models predic¬ 
tions in both cases are comparable with the computations of the 
comprehensive model. These results indicate that the new models 
introduced in this paper can be effectively used for practical design 
purposes. 

The validation shown in Figs. 9 and 10 reveals that the time of 
pyrolysis, an important parameter for reactor engineering, can be 
well predicted by the simplified models. So, it can be concluded that 
the formulation presented for the thermally thin particle provides 
a useful modeling tool which can be as accurate as the comprehen¬ 
sive model. The main advantage of the simplified models compared 
to the comprehensive model is that it can be easily implemented 
into a reactor model and solved much quicker. In industrial reac¬ 
tors such as fluidized beds and entrained-flow reactors, a designer 
needs to deal with thermochemical conversion of a large number 
of small particles. The most important parameters to be calculated 
from a single particle model are the history of temperature and net 
heat flux at the surface, and mass loss. The results of this study show 
that the time and space integral model can be treated as a useful 
tool for engineering applications. 

In real furnaces, the shape of particles could be far from a slab 
and they may be cylindrical and/or spherical-like. Therefore, it 
is worth extending the analysis of the time and space integral 
method to the shapes other than a slab. The derivation of the 
key relationships would be very much similar to that presented 
for a slab particle. However, one should realize that unlike slab 



Time [sec] 


Fig. 8. Comparison of the prediction of the simple models with the experiment and 
the prediction of Lu et al. model [4] (1 = 160 jjim and T r = 16251<). 


particles in that the heat transfer area is identical along the par¬ 
ticle, in a cylindrical and/or spherical particle the heat transfer 
area consistently decreases from its maximum values being at the 
surface to the center of particle where the surface area is zero. 
This has to be accounted for when performing space integration 



Half thickness (pm] 



Heat flux (kW/nr) 

Fig. 9. Conversion time of beech wood particle predicted by a comprehensive model 
and the simplified model: (a) conversion time versus particle half thickness q" ext = 
100 kW/m 2 ; (b) conversion time of a 1 mm particle versus heat flux. 
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Fig. 10. Conversion time of spruce wood particle predicted by a comprehensive 
model and the simplified model: (a) conversion time versus particle half thickness 
q" ext = 100 l<W/m 2 ; (b) conversion time of a 1 mm particle versus heat flux. 


of the heat transfer equations in all stages of the process described 
previously. 

Below, we present the expressions for computation of time of 
pyrolysis initiation, pyrolysis time and post-pyrolysis heating up 
duration suitable for cylindrical/spherical particles. The objective is 
to provide plant engineers and designers with simple relationships 
to help them estimate preheating and pyrolysis time as well as mass 
loss history. The form of the results for cylindrical and spherical 
particles is the same as those of slab, but with different numerical 
constants. In general, the relationships for evaluating the durations 
of preheating, pyrolysis and post-pyrolysis heating are as follows: 


fp.ini — 


R 2 


C 2 k B R(T p - T sR ) - (C 2 /C, )(q" - q")R 2 


£p — fp.ini + 


R 2 {PcCpCQp + C 2 PBCpBQp } 


p.ini' 


Cikc(q'p + C ini ) 


R 2 (C 2 T r -T sp C 2 
ph a c \ 2 T sp - T p + C, 


(43) 


(44) 


Cia B (l +((7 o/<?r)) 


«b(<?p + (Jr) 


where for a cylinder (Ci =4, C 2 = l) and a sphere (Ci =5, C 2 = ^/3). 
Looking back into Eqs. (15), (32) and (39) one may realize that for 
a slab Ci =3 and C 2 = 2. Notice that the relationship for estimating 
the pyrolysis time; i.e. Eq. (43), is obtained based on the single¬ 
temperature profile approach. 

The accuracy of this model for prediction of the mass loss of 
cylindrical and spherical particles is assessed using the data of 
Lu et al. [4], who investigated the effect of particle shape on the 
pyrolysis of small particles. In their study, cylindrical and spherical- 
like sawdust particles with equivalent diameter of 320 p,m and 
nearly the same mass and volume were pyrolyzed at 1625 K. The 
aspect ratios of the cylindrical-like and spherical-like particles were 
6 and 1.65, respectively. Fig. 11 compares the prediction of the 
simple model with the experiments and the model prediction of 
Lu et al. [4] for both cylindrical and near-spherical particles. As 
evident from Fig. 11, in both cases the simple model prediction 
compares well quantitatively and qualitatively with the prediction 
of a comprehensive model and the measured mass loss history. A 
further observation from this figure is that the spherical-like par¬ 
ticle pyrolyzes slower than the cylindrical particle does, which is 
due to the lower surface area-to-volume of the near-spherical par¬ 
ticle than the cylindrical sawdust species. Lu et al. [4] developed a 
3-D particle shape reconstruction algorithm to measure and cal¬ 
culate particle surface area and volume. The measured external 
surface area for slab-like, cylindrical and near-spherical particles 
was 0.491 mm 2 , 0.479 mm 2 and 0.344 mm 2 , respectively. 

We place an emphasis on the key parameter which needs to 
be carefully assigned a constant value; i.e. the pyrolysis tempera¬ 
ture T p . Despite that this parameter may not exactly represent the 
real temperature at which the decomposition takes place, it greatly 
facilitates simplification of a pyrolysis model compared to the case 
in which it is calculated from a finite-rate kinetics; e.g. Ref. [25,26]. 
Given that a wide range of values has been reported in the liter¬ 
ature, T p should be treated as an adjustable parameter leading to 
results comparable with experiments or prediction of a detailed 
model based on finite rate kinetics [28]. 

Past studies by Spearpoint and Quintiere [23,32] indicated that 
T p obtained from experimental tests is material and orientation 
dependent. They discovered that when the external heat flux was 
along the grain, the ignition temperatures were much greater than 
those obtained for the case that the incident heat flux was radiated 
across the grain. For instance, they determined average ignition 
temperatures of 375 °C and 204 °C for Redwood for the external 
heat flux radiating along and across the grain, respectively. Another 
study by Moghtaderi et al. [33] revealed that the ignition tempera¬ 
ture depends on the external heat flux and moisture content. They 
measured ignition temperatures for Radiata pine wood for inci¬ 
dent heat flux of 20-60 kW/m 2 and moisture content of 0-30%, and 
found a range of 298 °C-400 °C. Moreover, Yang et al. [34] reported 
a range of 190°C-310°C for the ignition temperature of wood. As 
a result, a designer needs to carefully select a proper value for the 
pyrolysis temperature otherwise it may lead to insufficient accu¬ 
racy of the simplified models predictions. If experimental data is 
available, it is suggested to assign a measured value for T p , other¬ 
wise the correlation proposed by Park et al. [35] may be employed 
as a useful tool which allows one to estimate T p as a function of 


(42) process parameters such external heat flux and particle size. 
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Fig. 11. Comparison of the prediction of the simple model with the experiment and 
the prediction of Lu et al. model [4]: (a) cylindrical-like particle; (b) spherical-like 
particle (d eq = 320 |jim, T r = 16251<). 


5. Conclusion 

A simple pyrolysis model for a thermally thin biomass particle 
is developed based on a time and space integral method, which has 
not been previously presented in the literature. The model is based 
on the main assumptions of constant thermo-physical properties, 
decomposition of biomass according to a shrinking core model at 
a thin reaction layer at a prescribed pyrolysis temperature, and 
approximating the spatial temperature profile with a polynomial 
function. The formulation of various stages of the conversion pro¬ 
cess has been presented using the time and space integral method 
resulting in transformation of partial differential form of heat trans¬ 
fer equation into an algebraic equation. 

Two different treatments were presented. The first formulation 
allows one to compute the history of key process parameters such 
as net heat flux at the surface, mass loss rate, char penetration depth 
and particle weight loss according to the calculation algorithm 
presented in Fig. 7. This model may, for example, be used in com¬ 
bustor and gasifier design codes where a large number of biomass 


particles undergo a thermal decomposition, since computation¬ 
ally it is cheaper and easier to implement in a CFD code than 
comprehensive models. On the other hand, the second treatment 
provides rather simple relationships for estimating the duration of 
various stages of the process including preheating, pyrolysis and 
post-pyrolysis heating. This method can be the interest of plant 
engineers since it provides a simple but useful tool to sufficiently 
predict the mass loss history of a pyrolyzing particle (see Fig. 5). 
Also, time evolution of the surface temperature can be approxi¬ 
mated using the method illustrated in Fig. 6. 

Typical experimental data related to the mass loss history of 
sawdust particles reported in the literature are used to examine 
the accuracy of the presented models. The predictions are found 
to be in a good agreement both quantitatively and qualitatively 
with the experimental data as well as the prediction of comprehen¬ 
sive numerical models. Additional validation is carried out using 
the prediction of the comprehensive pyrolysis model previously 
published by the authors [1,2] in terms of the conversion time of 
beech and spruce wood particles. The predictions of both treat¬ 
ments based on the time and space integral model agree well with 
the results of the comprehensive model. The key conclusion is that 
if the thermo-physical properties are assigned proper values, the 
pyrolysis models based on double- and single-temperature profile 
can be effectively used in practical applications. 
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